Non-uniform vortex lattices in inhomogeneous rotating Bose-Einstein condensates 
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We derive a general framework, in terms of elastic theory, for describing the distortion of the 
vortex lattice in a rotating Bose-Einstein condensate at arbitrary rotation speed and determining 
the dependence of the distortion on the density inhomogeneity of the system. In the rapidly rotating 
limit, we derive the energetics in terms of Landau levels, including excitation to higher levels; the 
distortion depends on the excitation of higher levels as well as on the density gradient. As we 
show, the dominant effect of higher Landau levels in a distorted lattice in equilibrium is simply to 
renormalize the frequency entering the lowest Landau level condensate wave function - from the 
transverse trap frequency, to, to the rotational frequency, fi, of the system. Finally, we show how the 
equilibrium lattice distortion emerges from elastohydrodynamic theory for inhomogeneous systems. 

PACS numbers: 05.30.Jp, 03.75.Hh, 03.75.Kk 



i 

-a 
c 

o 
o 



> 

(N 
00 

o 

\o 
o 

-I— > 



i 

C 

O 

o 



X 



I. INTRODUCTION 

Studies of vortices in rapidly rotating Bose-Einstein 
condensates [| have enabled one to look in detail at the 
structure of the vortex lattice. While, to a first approx- 
imation the vortices form a uniform triangular lattice in 
the limit of a large number of vortices, N v , Refs. and 
§ have shown, for slow rotation, that the vortex lat- 
tice undergoes a slight distortion, becoming more widely 
spaced near the edge of the cloud. The distortion ob- 
served later by Coddington et al. Q is in good agreement 
with theoretical predictions. The distortion expected in 
the fast rotation - lowest Landau level (LLL) - limit in 
a harmonic trap was discussed analytically in Ref. |5j 
and numerically in Ref. [(|. Our aim in this paper is to 
produce a general framework for calculating vortex lat- 
tice distortions at arbitrary rotation rates in terms of the 
energetics of the distorted lattice. 

The distortion of the vortex lattice is driven by the 
particle-density gradient in the system. Since the energy 
of a vortex increases with increasing local density, vor- 
tices tend to be forced towards regions of lower density. 
For clarity, we focus on the directions transverse to the 
rotation axis (z), and assume the system to be uniform 
along the rotation axis; the ready generalization to in- 
clude structure in the direction of the rotation axis does 
not change the results conceptually. In Sec. II we give a 
general discussion of the energy of a non-uniform vortex 
lattice. This is applied to the case of slow rotation in Sec. 
Ill and to the case of rotation rates comparable to the 
trap frequency in Sec. IV. Finally, in Sec. V we explore 
the relationship between the present work and elastohy- 
drodynamic theory 0, H, U E3 which describes the long 
wavelength modes of the vortex lattice. We take 1% = 1 
throughout. 



II. ENERGETICS OF VORTEX 
DISPLACEMENTS 

We assume that the vortex displacements from a uni- 
form lattice vary slowly in space and are described by a 
smooth displacement field, e(r), and that the lattice ro- 
tates at angular frequency f2. We work primarily in the 
rotating frame, where the total energy, E 1 = E — flL z , is 
generally a functional of the independent degrees of free- 
dom - the displacements of the vortices, the smoothed 
density, n(r), and the smoothed superfluid velocity in 
the rotating frame, which we denote by vr. [We denote 
the local density, with all its wiggles near the vortices, by 
rt(r), and the microscopic superfluid velocity by v.] In 
the final section we consider the energy as a functional of 
all these variables. However, in equilibrium, the super- 
fluid velocity is determined by the displacement of the 
vortices by the equation for quantization of circulation. 
In this case it is most convenient to regard the energy 
E', for a slow variation of the lattice distortion, as a 
functional only of the smoothed density, n(r), and the 
smoothed vortex density, n v (r). The important terms, 
as we shall see, are the first order change in the energy 
with e, which is present only for a spatially non- uniform 
system, and the change in the kinetic energy of rotation 
induced by a lattice distortion, which leads to a term 
second order in e. 

Since a uniform triangular lattice does not minimize 
E 1 ', the first order variation of E' under a distortion of 
the vortex lattice is non- vanishing and of the form 



5E' 



d 2 r Sn v (r)Q(n,n v ) 



(1) 



(The energy Q is the analog for a vortex of the quasipar- 
ticle energy in the Landau theory of Fermi liquids.) To 
first order in e, 



Sn v (r) = -n°V • e(r) , 



(2) 
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where n° = mSl/n is the vortex density of the uniform 
lattice. Integrating by parts, we have 



8^E' = 



VQ, 



(3) 



where we henceforth denote the integration f d r over 
the transverse directions simply by J, and to first order 
in e we replace the n v in Q by n® . (The second order con- 
tribution from the dependence of Q on n v leads to terms 
of relative order 1/N V [TJ.) Since VQ ~ (dQ/dn)Vn, 
and dQ/dn is generally positive, vortices are driven to- 
wards regions with lower density. They are prevented 
from moving too far, however, by the rotational kinetic 
energy of the system which increases as they move away 
from a perfect triangular lattice. 

To determine the kinetic energy of flow in terms of e, 
we start with the equation for quantization of circula- 
tion to relate the modification of the flow velocity to the 
vortex displacements: 



I v(r') ■ di 

/cm 



2tt 



N v (r), 



(4) 



where v is the microscopic local velocity in the lab frame, 
the contour C is a circle of radius r about the origin, 
and N v (r) is the number of vortices contained within 
the circle. For a uniform lattice with N v (r) replaced by 
irn®r 2 , the flow is solid body with v<f, = fir. Under a 
local radial displacement of the vortex lattice, SN v (r) = 
— 27re r rn° at fixed r. Thus the smooth flow of the fluid 
in the azimuthal direction is changed by changed by 



2tt 

Vr,4, = n v e r , 

m 



(5) 



which is the fluid flow in the rotating frame. An outward 
displacement of the vortices leads to a slowing down of 
the rotational flow. The kinetic energy of the smoothed 
flow in the rotating frame is 



K 



— I n(r)mv 2 — O j n(r)mrv^ 



n(r)m(v — rfl(py 



-in 2 



(G) 



where the integral in the second term of the first line is 
the angular momentum of the flow (the energy associated 
with the variations about smooth flow is included in Q); 
and I = J mn(r)r 2 is the moment of inertia calculated for 
the smoothed density profile. A distortion of the vortex 
lattice at fixed n(r), leads to a second order contribution 
to the kinetic energy 



SK = - j mn(r)v 2 R ^, 



(7) 



which, with JSJ, becomes the second order contribution 
to E' in terms of e: 



&>Ef = I 2 -^(n°) 2 n(r)e 2 r . 
J m 



(8) 



Corrections to this result are of relative order l/N v . 

The energy to second order also includes contributions 
arising from gradients of e, and this has the usual form 
for an elastic medium, 



2Ci(V 



Co 



dx 



dey 
dy 



de x 
dy 



d€y 

dx 



(9) 



where C\{n) is the compressional modulus and C2{n) the 
shear modulus of the vortex lattice 01 . Since the C's 
are of order fin, and Ve is of order e/R, where R is the 
radius of the system, the elastic energy density is of order 
nne 2 /R 2 . By contrast the kinetic energy density in (JSJ is 
of order mf2 2 ne 2 (where £ = l/^/irn v is the characteristic 
spacing of vortices), a factor ~ R 2 ji 2 — N v greater than 
the elastic energy ©; thus the elastic energy may be 
neglected in determining the distortion of the lattice. 

The two important terms in the change in the energy 
under a vortex displacement together are 



SE' = n°e- VQ(n) 



2n< 



(n° v ) 2 n(r)e 2 r . (10) 



Minimizing SE' with respect to e r we find the equilibrium 
lattice distortion fl4j |. 



1 dQ 

inn dr 



(11) 



As we now show, Q ~ n/m in the slow and fast rotation 
limits, so that e ~ ~£ 2 dlnn/dr. 



III. SLOW ROTATION 

In the slowly rotating limit, the kinetic energy of 
an individual vortex at position Rj is approximately 
J. n(r)/2mr 2 « 'yn(Rj), where n is the microscopic par- 
ticle density in the region of the vortex, the integra- 
tion is over the Wigner-Seitz cell of radius £ containing 
the vortex, r is measured from the center of the vortex, 
7 = (7r/m) ln(£/£), £ is proportional to the vortex core 
size, and the smoothed density, n, is evaluated at the po- 
sition of the vortex. Thus the leading term in the energy 
of the vortices is 



E mT t = J n v (r)n(r)j. (12) 

To leading logarithms, Q = rvy in this limit, and 

^ 2 , /„ dhxnir) , „. 

£r ^4 ln(W) * ' (13) 

in agreement with Refs. i and 0. 

Before turning to the rapid rotation limit, we remark 
that one can calculate Q in general between the slow ro- 
tation and very rapid rotation limits using the approach 
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of Ref. The calculations, which we do not do here, 
are complicated by the presence of explicit non-negligible 
contributions to the vortex energy that depend on dn/dr 
as well as d 2 n/dr 2 , as well as by the need to include 
explicitly the local fluid velocities in the intermediate 
regime. 



IV. LOWEST LANDAU LEVEL LIMIT 

A Bose-condensed system rotating at frequency close 
to the transverse trap frequency, uj, can be described by 
the trial condensate wave function, 



*(r) = N 1/2 h{r)ct> 



LLL, 



where 0lll = X e ~ r ^ 2d is composed only of lowest Lan- 
dau levels, with \ a polynomial in z = x + iy; d 2 = 1/mu), 
and N is the total number of particles. As long as the in- 
teraction energy per particle is small compared with huj, 
one does not expect the modification h of the LLL wave 
function, which admixes higher Landau levels, to vary 
significantly over the scale of the intervortex separation. 
It is thus a reasonable approximation to assume that h 
is a real slowly varying function dependent only on r. 
We first derive the angular momentum in the state 



N I d 2 r X *h 2 e- r2 / d2 z^, 
J oz 



(15) 



where in the latter form we use the fact that (zd/dz 
z*d / dz*)r 2 — 0. Then integrating by parts we have 



(L) = -N j ^L(z X *h 2 e-r 2 ' d2 )x 

d\nh(r) 



-N + uil — I fir 



dr 



(16) 



where n(r) = ^(r)! 2 and I — m J fir 2 is the moment 
of inertia. [The final term in Eq. Ijl6|l corrects the ex- 
pression for (L) in Ref. [Bj].] More generally, for real 
h(r) = h(x,y), the expectation value of the angular mo- 
mentum is 



( L ) = I M(^-l)-r.VlnMr) 



= -N + ujI- 



J h r • V 



ln/i(r). (17) 



To calculate the expectation value of the Hamiltonian 
in the state ljTl|l , 

E = I * ( + ^mcj 2 rA * + E int , (18) 



where E- luii is the energy due to interparticle interactions, 
we note that for a lowest Landau level, 



w2 i 1 2 2, 

d , d 



-r 2 /2d 2 
r 2 /2d 2 



(19) 



By a calculation similar to that of (L) above, we find [Uj . 



E = uj I + / - — 



h ( d In h 



2m \ dr 
dlnh(r) 
dr 
(14) 

For a more general h(r), the energy is 



—u> / nr 



E; 



int • 



(20) 



E = uj 2 I+ / — [VlnWr)] 2 - uj / ri r • Vlnft(r) 
J 2m J 

+E int .(21) 

Since dlnh(r)/dr is of order 1/i?, we may replace n(r) 
by its coarse grained average n(r) in the terms in E and 
(L) involving d\nh(r)/dr. The energy E' = E — fl(L) in 
the rotating frame is then 



E' = uj{uj - n)I + QN + 
— (uj — f2) / n(r)r 



n(r) /din h 



2m \ dr 
dlnh(r) bg I 2 



dr 



b 9 f 
— J n(r) 



(22) 



where b is the Abrikosov parameter. 

To determine the structure of the condensate it is most 
convenient to regard h(r) and the smoothed density n(r) 
as the independent variables in the energy (|22|l . Varia- 
tions of h(r) at fixed n(r) must change the local vortex 
density, n v (r). To see how, we write t he g eneralization 
for the state l|14(l of the result of Ref. [lj relating the 
particle density to the vortex positions: 



Vh/ n(r) 

4 



h 2 (r 



= —muj + nn v (r). 



Thus, under a variation of h(r) at fixed n(r), 

Sn v (r) = -— 5V 2 ln/i(r). 
2tt 



(23) 



(24) 



While the moment of inertia, /, for the smoothed density 
is independent of h(r) at fixed n(r), the difference I — I 
does dep end on the density of vortices, and, as shown in 
Ref. [Tl|, has the structure, 



I — I 



m n(r) 
7r n v (r) 



(25) 



Thus under a variation of h(r), we find, after integration 
by parts, that 



5(1-1) 



1 dn(r) d\nh 
-o- 



if2 2 dr dr 



(26) 



4 



This term is of relative order X/N v compared with the 
variation of the term — (u> — Q) J n(r)rdhi h(r)/dr in E' , 
and can be neglected. 

We now minimize the energy (12211 with respect to h(r) 
at fixed smooth density n(r) and find, 

Vlnh(r) = m(uj-Q)r, (27) 

so that 

h{r) = Ce m{ui - n)r2/2 . (28) 
Substituting this result back into Eq. (|14ll . we obtain 

*(r) ~ jv 1 / 2 xe~ mnr2/2 - (29) 

Although in terms of the frequency u>, this wave function 
includes higher Landau levels via h, their only effect is 
to change the oscillator frequency in the levels to the 
rotation frequency, £1. 

Minimization of E' with respect to n(r), with Eq. 1|28J) . 
yields the Thomas-Fermi profile, 



and thus 



1 



n{r) = Vg 2 



»(()) ( 1 



r 

~R- 



(30) 



plus terms of relative order 1 /N v . Here fi is the chemical 
potential. The radius of the cloud is R = [%Nbg / irni(u) 2 — 
fi 2 )] 1 / 4 . The Thomas-Fermi profile reflects the reduction 
of the effective trapping potential by the centrifugal po- 
tential. 

We now derive the distortion of the lattice. Measur- 
ing the displacements with respect to a uniform triangu- 
lar lattice of vortex density n^J = mfl/ir, and using use 
Eq. l|27 |) in ffi ty . with ©, we find, 

\ ^ ln = ~ m{U ~ n) ~ U ° V ' 6 (r) - (31) 

For equilibrium lattices, the displacement e is entirely in 
the radial direction, and we may integrate (|31|) to find, 



d\nh(r) 1 d\nn(r) 



+ m{ui — n)r + 2irn° v e r . (32) 



dr 2 dr 

With Eq. (|27J) for h, this result reduces to 

1 d ln n 



Amfl dr 



(33) 



which has the form with Q = irn/m. Using the 
Thomas-Fermi profile, we find 



1 



2mO R 2 — r 2 ' 



(34) 



Were we to assume a uniform lattice (e = 0), then 
Eq. would imply 



h(r) - ^e™'"- 1 ^ 2 / 2 , 



(35) 



(36) 



The admixture of higher Landau levels would not only 
modify the effective frequency, but modulate the wave 
function by the square root of the smoothed density. This 
result for h{r) does not minimize the energy, due to the 
imposed constraint of a uniform triangular lattice [l6| . 



A. Elastic energy in the LLL 

We now compute the terms in the energy associ- 
ated with the distortion of the vortex lattice. Inserting 
Eq. ^ for dh/dr into Eq. we derive the energy in 
the rotating frame, written in terms of the coarse grained 
density, n(r), and the lattice displacement: 



E'{n(r),e r (r)} = 



-I + lu{uj -0)(J - I) + Ei 



hit 



t r ^- + 2mn 2 J ne 2 r . (37) 

The penultimate term has the same form as that in the 
slowly rotating limit, (|12(l . only with 7 = ir/m. The fi- 
nal term is the same as in (|1C)|> . and is the modification 
of the kinetic energy associated with the decreased az- 
imuthal velocity caused by radial vortex displacements. 
Again, minimizing Eq. (|37J) with respect to e r , dropping 
the corrections arising from the difference between / and 
/, which are suppressed by a factor 1 — fl/tu, we are led 
directly to l|33|l . The angular momentum in terms of e r 
and n is 

(L) =uj(I-I) + m-2mVt [ nre r . (38) 



V. CONNECTION WITH 
ELASTOHYDRO DYNAMICS 

The results for the equilibrium lattice distortion are all 
contained in the elastohydrodynamic theory of Refs. 0, 
IE IE EH i when the dependence of the elastic energy on 
density gradients is taken into account. The key equa- 
tions in the elastohydrodynamic theory (all in the rotat- 
ing frame, denoted by a subscript R) are the conservation 
of circulation. 



v R + 2f2 x e = V$/m, 



(39) 



where $ is the superfluid phase; the superfluid accelera- 
tion equation, the time derivative of Eq. i|39|) , 



—— + 2f2 x e 

at 



-V(M + Voff), 



(40) 



where fj, is the chemical potential of the matter; and the 
equation for conservation of momentum, 



rn- 



dt 



2mfl xj R + VP + nWcff 



(41) 
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Here v# and jn are the smoothed superfluid velocity and 
current density in the rotating frame, P is the pressure, 
and V e g(r) = Vt rap (r) — \m£l 2 r 2 . The elastic force, — <x, 
is given by 



/ \ SE' 
* (r,t) = 

oe 



(42) 



where E 1 is the full energy in the rotating frame for the 
distorted lattice, including the usual elastic energy 
In the elastohydrodynamics, which also describes non- 
equilibrium dynamics, the smoothed density, n(r), the 
displacement, e(r) , and the smoothed superfluid veloc- 
ity, Vfl(r), are independent dynamical degrees of free- 
dom, and one must regard E' as a functional of these 
three independent variables. Only in the determination 
of e in equilibrium can one eliminate the superfluid veloc- 
ity as an independent variable. Keeping v# fixed means 
that the kinetic energy variation, 5K, and hence S^ 2 'E', 
Eqs. and (jHJ, do not contribute to cr. 

In equilibrium, where one has axial symmetry, Eq. i|39|) 
reduces simply to ©; then from Eq. ijiOjl. d(/j,+V e s) = 0. 
In addition, at zero temperature, VP = nV/i. 

To show how the displacement of the equilibrium lat- 
tice, i|S3j, emerges from the elastohydrodynamics, we 
write jij = nvR and subtract Eq. (|40|l from (|41|l divided 
by n, to find |l7j . 



The radial component of this equation in equilibrium 
reads 



(44) 



so that from Eq. JSJ), 



2mfi x (v R — e) = . 

n 



(43) 



4mfi 2 £ r = -. (45) 

n 



Evaluating a from the elastic energy contribution 
Q Je r dn/dr in Eq. (|37|l we are led immediately to the 
vortex displacement (|3"5|) in the LLL. 
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